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Abstract: We propose a simple model for studying small-x physics in which we take 
only the coUinearly enhanced part of leading and subleading kernels, for all possible 
transverse momentum orderings. The small-x equation reduces to a second order 
differential equation in t = logfc^/A^ space, whose perturbative and strong-coupling 
features are investigated both analytically and numerically. For two-scale processes, 
we clarify the transition mechanism between the perturbative, non Regge regime and 
the strong-coupling Pomeron behavior. 
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1. Introduction 



The problem of subleading logs corrections to the BFKL equation [1-4] has been re- 
cently investigated by several authors [5-11] in connection with the small-x behavior 
of structure functions and of two-scale processes [12-14]. 

We have advocated [8,9,11] in this context, the importance of incorporating 
Renormalisation Group (RG) improvements and collinear physics in a new form of 
the small-a; equation. The latter has the virtue of including all collinearly enhanced 
higher order contributions to the kernel, and provides stable estimates of the re- 
summed gluon anomalous dimension and of the hard Pomeron. 

Here we present a simple, but powerful tool for studying the problem of small-x 
physics that we call the collinear model — namely a model where all and only the 
coUinearly-enhanced physics is correctly included, in particular the full dependence 
on the one-loop running coupling, the splitting functions, and the so-called energy 
scale terms [4, 8]. 

The model has the properties of correctly reproducing the one-loop renormali- 
sation group results and of being symmetric (given a problem with two transverse 
scales, its results do not depend on which of the scales is larger). These are properties 
desired from the resummation of the NLL corrections in the case of the full BFKL 
kernel. 

While it does not correctly resum the series of leading and subleading logarithms 
of s (i.e. the non-coUinear part of the problem), in that region it does have a structure 
which qualitatively is very similar to that of the BFKL equation, and can usefully 
serve as a model. In contrast to the BFKL equation, it is very easily soluble, as a 
Schrodinger-like problem. 

2. The collinear model 

Let us start recalling [4] that, in order to specify a small-x model for a hard process 
with two scales k and ko, one should specify the energy scale so, the kernel and the 
impact factors. They all enter the fc-factorised form of the cross section 

(TABiY; K ko)= fP: ( hA{k) g^{k, k,) hB{ko) , (2.1) 

J 27ri \so{k,ko)J 

where Y = log s/sq and 

^, = [u;-/Cj-^ (2.2) 

is the gluon Green's function. In the following, we refer normally to the symmetrical 
factorised scale Sq — kkg, by keeping in mind that for k ':$> kg (or viceversa), one can 
switch to the relevant Bjorken variable k^/ s (or k^/s) by the similarity transformation 
induced by {k/koY on the kernel and the Green's function. 
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2.1 Definition of the model 



Our collinear model is then defined — in logarithmic variables t = log/c^/A^ — by 
specifying the kernel ICco{<ysit),t,t') whose «s(t) expansion corresponds to higher and 
higher order collinear singular kernel for both k ^ k' and k' ^ k. 





Figure 1: Subleading kernels being resummed by IC^ik^k') in Eq. (2.3). Wavy (dashed) 
lines denote high-energy (low-energy) gluon exchanges, corresponding to the Ois/u) {Aias) 
part of the gluon anomalous dimension. 

In other words, our model is such that it reproduces the DGLAP limits [15] 
for branchings with ordered transverse momenta, and the anti-DGLAP limit for 
branchings with anti-ordered transverse momenta. Accordingly, the kernel governing 
a small- a; branching of a gluon with transverse momentum k' into a gluon with 
transverse momentum k is (Fig. 1) 



1 hh' 

Lxut,t') = —icuk,k') 

LO LO 



(2.3) 



exp 



|-^(i - t') + A^ico) J^^ a,{r) drj e{t - t') + (t ^ t') . 



LO I 2 

The term proportional to Ai[lo) comes from the summation of the non-singular part 
of the DGLAP splitting function, including a full treatment of the running coupling: 

Nr.a, ^ l_ 
~ bt 



-fgg{as,Lo) = Ai{Lo)as{t) + 0{a^) , 

LO 



a. = 



TT 



(2.4) 



Therefore the kernel for a single small-x branching actually resums many branchings, 
of which the last (and only the last) is governed by the l/co part of the splitting 
function. 

The scale of «s is the larger of the two scales involved in the branching, as implied 
by DGLAP. Finally, the factor kk' is due to the fact that we work in t-space instead 
of fc-space, and the remaining one exp[—^Lj{t — t')] is due to the choice of scaling 
variable kko/s, as explained before. 

The collinear properties of the kernel (2.3) can also be seen in 7-space, by the 
expansion 



ic^Kk') 



J2[^,{e)r'K{k,k'), 

n=0 



(2.5) 
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where the scale- invariant kernels have eigenfunctions e^''' 2)* and eigenvalue func- 
tions Xn(7) given by [9] 



_ l-A,{A, + b)---{A, + {n-l)b) 1■{A^- b){A, - 26) • • • {A^ - nb) 

Xn W j / . 1 s „ j 1 ~r 



(7 + |a;)"+i 



(l-7+|a;)"+i 



/ 



Of particular interest is the NL expansion 

~ q;s(^) 



(2.6) 



(2.7) 



^0 - o'^ 1 :^ 



(1-7)^ 



+ as(i)M + NNL 



whose ci;-dependence can be reabsorbed [11] in a redefinition of impact factors and 
kernel, to yield the leading and next-to-leading BFKL kernels of the model: 



1 1 

Xo = - + 

Ai ^ A^-b Xo f I 
Xi = — 
T 



(2.8a) 
(2.8b) 



:i-7)2 2 V7^ (1-7)^ 
where the renormalisation scale for ctg has been taken as t, i.e. 

3CBFKL(as(^),7) = «s(^)Xo(7) + asXi(7) ■ 



At leading order, this reproduces the poles at 7 = and 1 of the true BFKL kernel. 
At next-to-leading order it reproduces the 7 = 0, 1 quadratic and cubic poles of the 
true BFKL kernel. Numerically the leading order kernel differs from BFKL quite 
significantly numerically, but retains a very similar structure — a saddle point at 
7 = 1/2, implying a power growth of the cross section, and diffusion. The double 
and triple-polar part of the ncxt-to-leading kernel turns out to be very close, even 
numerically, to the full BFKL NLO kernel, reproducing it to better than 7% over 
the whole range of 7 from to 1 (note though that this coUincar kernel has single 
polar and other parts, so that as a whole it may not be quite this close to the full 
BFKL NLO kernel). This suggests that coUinearly enhanced effects dominate the 
NLO kernel. 



2.2 First order formulation 

The main advantage of our coUinear kernel from the point of view of this article is 
its relative simplicity. Specifically it can be written in factorised form: 

X^{t, t') = U{t) V{t') Q{t - t') + U{t') V{t) Q{t' - t) , (2.9) 
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where 



U{t) ^a,{t)e^p^-^-^t + Ai{u) J as{r)dr^ , (2.10a) 
=exp|i^t- AH y"*a,(T)dT| . (2.10b) 
This allows us to recast the homogeneous BFKL equation^, 

/t POO 
dt' V{t')3'{t') + V{t) / dt' U{t')3'{t') (2.11) 
-oo Jt 

as a differential equation. Dividing IT into two parts, 

A{t)^U{t) I dt' V{t')'J{t') (2.12a) 

J — oo 

/oo 
dt' U{t') S'it') (2.12b) 

ouS'it) = Ait) + B{t) (2.12c) 

and taking tlie derivative leads to a pair of coupled differential equations: 

dA W 

'^ = ^A + UV3^, (2.13a) 
dB V 

— ^—B-UV^. (2.13b) 
dt V ^ ' 

For the specific kernel (2.3) that we consider, we have 

d^ f l+uj . oi'\ , _ . . X 

— = — + Aias + ^M + «s3^ 2.14a 

dt V 2 as J 

f = - Ai«s) B - a. J (2.14b) 

Since we have two coupled equations, there are two independent solutions. Examin- 
ing the equation for large and positive t, where ag is small, one sees that they can 
be classified as a regular solution 

:rK~exp(^-^-^t^, (2.15) 

which is dominated by A, and an irregular solution 



3^1 - exp ( ^-^t ] , (2.16) 



dominated by B. 



^ Since we work in f-space, the density 3^{t) differs by a factor of k from the customary [11] BFKL 
solution J^{k), at scale kko- 
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2.3 Second order formulation 

The coupled set of differential equations (2.13) can be recast in the form of a simple 
second order equation for 9^. In fact, by using (2.12), we can first rewrite (2.13) in 
the form 



UVS". (2.17) 



Then, in order to eliminate the resolvents appearing in (2.17) we introduce the op- 
erator 

and the wronskian 

w{t) = W[U, V] = UV - U'V = a,{t) (^l+LU- 2Aia, - . (2.19) 

By applying to (2.17) we finally obtain 

cu-DtS^ =(---] UV3^ = -wm, (2.20) 

which is the second order formulation that we were looking for. 

Equation (2.20) can be recast in normal form by the similarity transformation 

3^ = const • ^/w{t) h{t) , (2.21) 

which leads to a Schrodinger type equation, 

{-dl + Kff) h = (2.22) 

I fw'V I fw'\ w U"V'-V"U' 
Kfr=T - - -- + ■ 2.23 

Note that the above derivation is valid for any form of the running coupling as{t) 
which extrapolates the perturbative form {bt)~^ into the strong coupling region 
around the Landau pole t = 0. 

In what follows we consider various regularisations of the Landau pole, in par- 
ticular: 

(a) aJt) ^ ^Q(t -t), (cutoff case), (2.24) 

bt 

(b) as{t) = ^e{t -t) + ^e{t-t), (frozen case), (2.25) 
where t > sets the boundary of the perturbative behavior. 



5 




Figure 2: Qualitative t-dependence of the effective potential for the regularisations of 
type (a) and (b) of the coupling strength. 

It is obvious that such different forms will change the form of the potential 
(Fig. 2) and thus the boundary conditions on 5" (or h) coming from the strong cou- 
pling region. A similar Schrodinger formulation was found [16] for the (Airy) diffusion 
model [17,18] with running coupling, with a potential which roughly corresponds to 
the bottom of the well in Fig. 2. 

We consider in particular the solution 3^R{t) (^'^(t)) of the homogeneous equation 
(2.20) which is regular for t +oo (t —> — oo). If both conditions are satisfied, 
3^R = 3^L is an eigenfunction of the BFKL equation. The pomeron singularity uj = ujp 
is the maximum value of cu for which this occurs (ground state) . 

li cu > u!]p, 3^R ^ and the two solutions have rather different features. Due to 
the locality of the differential equation, is independent of the regularisation in the 
region t > t. On the other hand, will be dependent on the behavior of as{t) for 
t < i. For instance, in the case of as{t) being frozen (2.25), has the exponential 
behavior 



3^L{t) e"', {t<t), 
where k — y^Ves(t < t) is found from (2.23) to be 



(2.26) 



-{1 + u;) - Aia,{i) 



-{l + u;)-A,a,{t) ^ 

Z UJ 



(2.27) 



On the other hand, if as{t) is cutoff, it is simpler to use directly the first-order 
formulation, (2.13). We note that we can write 



uj-Jl = Al + Bl = a,{t)aL{t) + BL{t) , 



(2.28) 
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where ai and are continuous aX t = t, and are exponentially behaved, ~ Bl 
exp(^(l + u!)t), for t <i. By using (2.14) we then obtain 



Ar 



which defines the boundary condition for the Al and Bl components of in this 
case. 

2.4 Factorisation of non-perturbative effects 

The basic tool for describing BFKL evolution is the Green's function Suj{t,to) = 
kkoQoj{k, ko), which satisfies the inhomogeneous small-x equation 

coSUt, to) = 5{t - to) + 9Ut: to), (2.30) 

and is supposed to be well-behaved for t, to ±oo. The problem of factorisation is 
the question of the (in) dependence on the non-perturbative strong-coupling region. 

For t to, 9ui satisfies the same differential equation as H and is thus a su- 
perposition of two independent solutions. The large-i behavior implies a regularity 
condition and suggests the expression 

%it, to) = 3'RmLito)Oit - to) + 3'LmRito)eito -t), {t^ to) (2.31) 

where {S^l) is the regular solution for t —>■ +oo {t —>■ — oo) defined in the previous 
subsection. 

Actually (2.31) is a rigorous consequence of the second-order formulation. In 
fact, Su) satisfies the differential equation 

ht + -) Su, = -DtS{t - to) . (2.32) 

With a little thought, one can realize that 9u; must contain a delta function term in 
the form 

%{t, to) - -8{t - to) + %{t, to) (2.33) 

UJ 

and hence 

ht + -)% = --^5{t-to) (2.34) 
\ UJ / uj-^ 

showing that S is continuous function t — to with discontinuous derivative. 
Eq. (2.34) is an inhomogeneous Schrodinger type equation with a delta source, and 
its solution can be found in standard textbooks to be just hke the RHS of Eq. (2.31) 
also for t — to- We conclude that 

S(t, to) = -S{t - to) + ?K(t)?L(to)e(t - to) + Mt)3'R{to)Q{to - t) (2.35) 

UJ 
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with the normahsation 



hR,L{t) 



W[hR,hL]=hRh'i^-h'j,hL^l. 



(2.36) 



U! 



The main consequence of (2.31) is that the regularisation dependence is factorised 
away in "Jl, whenever t or to are large enough. This happens in particular in the 
coUinear limit t — to ^ 1 relevant for structure functions. 

3. Solutions: analytical features 

The coUinear model just defined can be solved in principle as a Schrodinger prob- 
lem by known analytical and numerical techniques and for both types of solutions 
occurring in the Green's function (2.2) (i.e., the left-regular and the right-regular 
ones) . 

The regular solution 5"^ is, for large t, perturbative, i.e., independent of the po- 
tential in the strong coupling region t < i, while the left-regular one is dependent 
on the strong coupling boundary conditions through a reflection coefficient >S'(a;) of 
the (S-matrix, which occurs in its expression for t > t: 



where 3^r (5"/) denote the regular (irregular) solution for t > t, with the normalisation 



This allows to rewrite the Green's function for t > to and a; > a;p in the form 



If both t,to » 1, but t — to — C(l), Eq. (3.2) is dominated by the first term, 
the second being suppressed exponentially in to- This term is, on the other hand, 
defined by boundary conditions for t, to — > +oo only, and is therefore independent 
of the strong coupling region. Its analytical and numerical form will be discussed in 
more detail in the following. 

The second term in Eq. (3.2) carries the regularisation dependence and contains 
the leading a;-singularities, in particular cup (Ref. [16]). In y = log(s/A;A;o) space, the 
sum in Eq. (3.2) defines two asymptotic regimes, as we shall see. 

3.1 Perturbative regime: a;-expansion and WKB limit 

Approximate solutions in the largc-t region can be found, as in the full small-x equa- 
tion, by the method of the 7-rcpresentation and a;-expansion. The regular solution 
is approximated by the expression: 



9^L{t)^3'i{t) + S{u;)3^R{t), 



(3.1) 



W[hR,hl]^l. 



(3.2) 




(3.3) 
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with 



d^X"" = xil, uj) = Xo (7) + <^ 



xr(7) ' 



where, for the colhnear model, 

xr(7) = 



1 



7 + 



1-7 



5^' 



+ 



Ai -6 



(3.4) 

(3.5a) 
(3.5b) 



(7 + ia;)2 (l-7+|a;)2- 

In [11] this representation was extensively studied, and shown to be a solution of 
the problem up to next-to-leading order and to all orders for the colhnear structure. 
But this was only guaranteed to work true for small values of cu (whereas for the 
continuation with the DGLAP anomalous dimensions it is useful to be able to access 
high u as well). Also there was no way of determining the coefficient of any higher- 
order error introduced by the procedure. A comparison of this representation with 
the exact solution, as is possible in the colhnear model, is therefore important (cf. 
Sec. 4.1). 

The expressions (3.3-3.5) can be further specialized in the large-t limit, where 
(3.3) is dominated by a saddle point at 7 = 7: 



bujt = x(7,^) ■ 



(3.6) 



The latter is related to the WKB approximation for solving the differential equation 
(2.22). In fact, on the basis of Eqs. (2.22) and (2.36), one can prove the asymptotic 
expansion 



2K{t) 



exp \ — I k{t) dr 



X 



1 + 



(3.7) 



where K,{t) is defined in terms of the effective potential (2.23) as 



K''{t) = V,fi{t)^-{l+u; 



2Ai 



bt 



2Ai 



bt 



2 
but 



+ 0( -r 



1 

(3.8) 



and is related to the saddle point value 7(0;, t) of Eq. (3.6) by 



7 = - ( 1 + J 



(3.9) 



It is interesting to note that in the colhnear model, because of the differential 
equation (2.22), the present method yields the irregular solution also, which is ob- 
tained by just changing the sign of the WKB momentum K,{t), i.e.. 



Mt)^-' 



2K{t) 



exp 



k{t) dr 



(3.10) 
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This is useful for evaluating the Green's function which according to Eqs. (3.2), (3.7) 
and (3.10) takes the form 



^WKB 



it, to) 



w{t)w{to) \^^(_ r (^)dr^ (3.11) 



2u;^\ K{t)K{to) 



to 

to 

{t>to), 



+S{u) exp 1^- (^J^\ J^^ k{t) drj 



where tg is the zero of k and S{u!) has to be determined from the boundary conditions 
on Sit t — t. 

3.2 Strong-coupling features 

For intermediate, small, and negative t's, the collinear model enters a second regime, 
where both kinds of solutions oscillate. Given the definitions of the effective potential 
in (2.23), the regime's boundary is roughly given by t < Xm/buj, where K,{t) becomes 
negative, and Xm is the minimum (in 7) of the effective eigenvalue function x(7,a;) 
defined in Eq. (3.4). 

A basic question concerning this regime is the spectrum of /C^^, which provides 
the (^-singularities of the Green's function, which in turn determine the high-energy 
behavior of the cross section. 

The leading ct;-singularity is the Pomeron cjp, i.e., the maximum uj value for which 
"Jl and "Jfi match each other, providing a zero energy bound state in the potential 
(2.23). The Pomeron properties are dependent on the strong coupling boundary 
conditions which are basically the frozen CKg and cut-off cases considered in Sec. 2. 

The first case is fairly simple, because of the fixed value of as for t < t. In this 
region, because of the replacement 



exp 



the solution is a plane wave whose momentum k is given in Eq. (2.27). Therefore, 
cup defines the boundary of the continuum spectrum, given by k — 0, or 

a;p = . (3.12) 

This is to be compared with the saddle point definition of the hard Pomeron [11] 
u!s{t), which is obtained by minimising the solution (as a function of 7) of 

c^. = «s(^Xol7)+^.^) ■ (3.13) 
Taking the b — form for Xi the minimum is at 7 = 1/2, and cUg is given by 

a;, = - r- (3.14) 
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which is identical to the true uj^. 

In [11] a second critical a;-exponent was proposed, termed uic-, corresponding to 
the position of the rightmost singularity of the anomalous dimension of the integrated 
gluon distribution. In the coUinear model the integrated gluon distribution can be 
defined as 



g^{t) = [ dt' exp \-^-^{t - t') +Ai [ dr a,{r) 

J-oo L ^ Jt> . 



m-^.m, (3.15) 



because the corresponding density (at energy-scale A;^) e^^^g^{t) can be shown to 
satisfy, in the collinear limit, the usual DGLAP equation with anomalous dimension 
+ Ai)as. The singularity of the effective anomalous dimension 

7 = ^ln^c.W (3.16) 

is at the point uj — uJc{t) where gu}{t) goes to zero, i.e., where A{t) goes to zero. 

This critical exponent ujc is related to uj-p for the cutoff case, in which as{t) — 

for t < t. In fact, since A{t) ~ as{t), we have A — ioi t < t; but we noticed 
A(t) 

in Sec. 2 that _ ~ Qwi't) is instead continuous at t — t, so that for t ^ t+, u-p 
satisfies the boundary condition 

AiU) = ^s(^"+) , y^O. (3.17) 

This means that u!c{t) < u>p, because the depth of the well, determined by cu < u!s{t), 
is to be further decreased below up in order to have A{t+) = 0. 

The relationships just found {uj^^'^™^^ = ujs{t), u^'^"^ > ujc{i)) represent two 
extreme cases of the boundary condition dependence of cup. If the strength as{t) > 
is positive but has intermediate size and shape for t <i, we expect in general that 

oucit) < a;p < (Jsit) , (3.18) 

i.e., the lower and upper bounds mentioned in Ref. [11]. 

3.3 High energy behavior and diffusion 

It is widely believed that a two-scale process — described by a small-x equation of 
BFKL type — is perturbative for large enough t and to, while it becomes a strong 
couphng process if the energy is so large as to allow diffusion to small values of t ~ 
(A;2 ~ A2) [19,20]. 

In the collinear model the Green's function has the explicit expression (3.2), in 
which the strong-coupling information is clearly embodied in the "^'-matrix coeffi- 
cient" S{u;). Therefore it allows a direct study of the relative importance of the 
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"perturbative" part 3^^ 5"/ and of the "strong-coupling" part 5"^, induced 
by diffusion through the boundary conditions ai t = t. 

For large t and to, takes the approximate WKB form (3.11), that we study 
in the special case t — to, so that 



S(Y; t,t)^j^^ [l + S(.) exp (-2 d.) 



(3.19) 



provided k^^{t) dr ^ 1. 

We now notice that if Y = \og{s/kko) = logs/A^ — t is not too large, the 
expression (3.19) is dominated by a saddle point at a; = a;, such that 



(3.20) 



2 {but - Xm) ' 



and therefore 



^(y,t)...(*) + ^ = f + x™ = 4(l+..-^)* . (3.21) 

so that Cd{Y, t) is not much different from uJs{t), the saddle point exponent mentioned 
before (Eq. (3.13)). Furthermore, at this saddle point the phase function takes the 
value {x'Iri - Xm/2) 



provided 



2, r 1 1 f ht\^ 



bt ~ u;;^ < 2y ( ^ y a;, ^ ~ as ^ . (3.23) 



In other words, the saddle point (3.21) is self consistent, i.e., it exists in the 
WKB region of Eq. (3.23), provided the effective parameter as^^Y <^ 1. 

The actual evaluation of (3.19) is now performed by distorting the uj contour as 
in Fig. 3 and picking up the saddle-point value of the background integral and the 
a;-poles: 

g(y; t, t) ~ , ^ e'^^W^ X [1 + 0{a^,Y^)] + e"^^ Rpe-^^+'^^' . (3.24) 

Here we have kept the leading Pomeron pole, with residue Rp, and neglected the 
nearby poles > a) = Ci;^ -f- , which are anyhow exponentially suppressed because 
of the large phase function O {{a^Y^)-^^'^) of Eq. (3.22)^. 

^Of course, the background integral contribution is equivalent to the one of the poles a;„ < w, 
as it can be seen by displacing the contour further to the left. 
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CO. 



^ ^0<X X X 



Figure 3: The contour of integration for the Green's function Su>iY;t,t). 



The corrections to Eq. (3.19), depending on the (small) parameter a^Y"^, can be 
explicitly evaluated from higher order WKB terms. More simply, we can use the fact 
that the effective u values are close to the minimum of %, in order to use the Airy 
form [16] of the Green's function, which in the cutoff case is given by the expression 



where the variables 



V Xm 



mo) - W^mo) 



f2buY fx 



— -t) 



(3.25) 



(3.26) 



evaluated at ^ cj, turn out to be both large parameters ~ t"^^^ of opposite sign. 
By using the expansion [21] 



Ai(OBi(0 
and the saddle-point value 



1 / ^ 



5 

32 



(3.27) 



, , . 2 t f bt 

we see that the correction is indeed of order a^Y^. The quantitative evaluation 
requires a careful treatment of a;-ffuctuations around u> — cUs = 1/2Y. By using 
Eq. (3.27) and the integral 

2iv^^ ' 15' 
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we finally obtain, for the perturbative part of the Green's function (3.24) the sub- 
asymptotic correction factor 

i + Y^^^lx'^HY'. (3.28) 

The latter turns out to coincide with the "non-Regge correction" to the "Regge 
exponent" uJs{t)Y found by other authors [7] in different but related contexts. 

We notice, however, that the true Regge contribution is the second term in 

Eq. (3.24), which is of strong-couphng type, with a t-independent and eventually 
leading exponent cup. The perturbative part, which dominates in the large-i limit, 
comes from the background integral and has no reason to be Regge behaved. 

Thus, the appearance of the parameter a^Y^ is Eq. (3.27) signals just the exis- 
tence of a "quantum" wavelength in the solution for u < ujs{t). When a^Y^ ~ 1 the 
saddle point breaks down and the solutions enter the small-t regime. 

It is non trivial, however, that in the well defined intermediate regime (3.23), the 
exponent uJs{t) with the corrections (3.28), appears to be an observable quantity. On 
the other hand, the exponent oodt) — the formal anomalous dimension singularity 
— does not directly appear in the Y dependence, because the oscillating behavior of 
the S^'s is masked by the onset of a;p dominance for 

Y>Yt= ^-i. (3.29) 

We conclude from Eqs. (3.19) and (3.24) that the two-scale Green's function 
shows a perturbative (non-Regge) regime where the exponent uJs{t)Y shows up with 
calculable corrections (Eq. (3.28)), provided the parameter afF"^ is small. Even 
before the latter gets large, at Y > Yt the Pomeron-dominated regime takes over, 
characterized by the regular solution, which is confined to the strong-coupling region 
of small fs. 



4. Numerical results 

Here we concentrate on a couple of numerical aspects. Firstly, the direct calculation 
in t-space of the regular solution: this allows to test the 7-representation method 
used for the full small-a; equation [11]. Secondly, the direct evolution in Y space 
of the Green's function, in a simplified case where the uj and Ai dependences are 
removed: this allows us to elucidate the transition region between the two asymptotic 
regimes just mentioned. 

4.1 Reguleir solution 

General structure. The regular solution J^r of the differential equation is ob- 
tained by starting at large t with an arbitrary initial condition and then evolving 
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downwards. Since the irregular solution falls rapidly relative to the regular solution 
when reducing t, one quickly reaches a situation where only the regular solution is 
left. By starting at large enough t one can ensure that this be true to arbitrary 
accuracy. 

In [11] much use was made of the (^-representation to obtain an approximation 
to the regular solution. In this section we wish to carry out several tests of the 
practical accuracy of the (^-representation. All results shown here were obtained 
using Ai = —1 and 6 = 1. 

First we illustrate the form of the unintegrated distribution in Fig. 4, normalised 
to be 1 at the maximum. Results are shown both from a direct solution of the 
differential equation and from the cu-representation. One sees that for larger values 
of t, the (^-representation is in good agreement with the exact solution, while for 
smaller t, where the solutions oscillate the results from the (^-representation are 
slightly out of phase with the exact solution. In general we are interested in the 
behavior to the right of the rightmost zero. 

1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 

Figure 4: Regular (unintegrated) solution from the explicit solution of the differential 
equation and from the u representation; shown for uj = 0.15. 

Critical exponent. One quantitative test of the ci;-representation concerns the 
critical exponent uJc, i.e. the value of u at which the anomalous dimension diverges. 
We recall that this is connected with the position of the rightmost zero of the regular 
solution: 3^u)c(^) = 0- The ci;- representation determination of the position of this 
zero, or equivalently its determination of uJc as a function of t, involves a small error 
which we call ScUc- Fig. 5 shows ScUc/cUc as a function of cuc, and we see that the 
relative error on cUc goes roughly as cu'^, or equivalently as a^. This corresponds to 
a NNL correction and is beyond our level of approximation. We note also that even 
for relatively large values oi cu ~ 0.3, the relative correction remains of the order of 



differential equation 

(o-representation 
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5% which is quite acceptable. In other words the NNL correction that arises is not 
accompanied by a large coefficient. 




Anomalous dimensions. Our second quantitative test of the CL^-representation 
concerns the anomalous dimension. The error in the a;-representation anomalous 
dimension, (^7, is plotted in Fig. 6 as a function of ctg for two values of lo. Let us first 
concentrate on the region for ctg < 0.01. We see that the error is roughly independent 
of a;, and proportional to o!^. In other words the difference between the exact result 
and the ^-representation is a term of O (fSg). We recall that the terms that we wish 
to include properly are the leading terms (oL^jujY' . the NL terms aj^a^jujy^ and the 
collinear terms {Q.^juS)uj^'' . A first correction of O {a^ is consistent with all these 
terms having been correctly included. 

For cis > 0.01 wc sec that the error in the anomalous dimension has a more 
complicated behavior: it changes sign (the dip) and then diverges, at the same point 
as the divergence in the anomalous dimension itself (solid curve): this is just a 
consequence of the position of the divergence of the anomalous dimension from the 
(^-representation being slightly different from the true position (cf. Fig. 5). 

4.2 Green's function with running coupling, diffusion, the infra-red 

We have already discussed in Sec. 3.3 the analytical expression of the Green's function 
of our model, and the features of the perturbative and strong-coupling regimes in the 
y-dependence. Here we want to study the transition region between the two, which 
is analytically hard to describe. 
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Figure 6: Error in the anomalous dimension from the a;-representation. 



Note first that several aspects of the regimes are due to just running ccg effects. 
The latter enters at two levels: first, because of the "acceleration" in the well of 
Fig. 2, it modifies the traditional fixed-Og kind of diffusion, weighting it towards 
lower transverse scales — this may well be connected with the breakdown of the 
saddle-point approximation for a^^'^Y ~ 1 as discussed above. Secondly, because of 
the strong-coupling boundary condition, it also introduces a qualitatively new kind of 
diffusion, perhaps more properly referred to as 'tunneling': namely there is a certain 
Y — Yt, defined in Eq. (3.29), at which the non-perturbative pomeron suddenly takes 
over, because 



and beyond which the Green's function is dominated by the regular solution, and 
thus confined to small i- values. 

This 'tunneling' phenomenon is qualitatively different from diffusion, in so far 
as there is not a gradual decrease in the relevant scale for the evolution, until non- 
perturbative scales become important, but rather there is a point beyond which low 
scales suddenly dominate, without intermediate scales ever having been relevant. To 
see this consider that the contribution to the evolution at scale t from an intermediate 
scale ti <^t (with a corresponding exponent uji ~ Xm/bti), is of order 

^XmY/{bti)-{t-ti) _ ^4 2) 

Prom this one sees, because of Eq. (4.1), that ti — t becomes relevant before higher 
scales do. 
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To study these effects at a qualitative level it suffices to consider a very simplified 
version of the coUinear model: one which retains only the running of the coupling, but 
not the a;-shifts of the 7 = 0, 1 poles, nor the Ai component of the NL corrections. 
We then examine the solution to 

g(y, t, to) = S{t - to)S{Y) + [ dyK® g(y) , (4.3) 

for this simplified kernel and we consider the effective exponent of the evolution 

c^efr = ^lnS(F,t,t) (4.4) 

as a function of Y . Fig. 7 illustrates the basic behavior of uj^fi for rather extreme 
kinematics — not intended to be phenomenologically relevant, but rather to show 
clearly the relevant features. Two values of the infra-red cutoff are considered. What 
is seen is that the exponent at first increases slowly and smoothly, and then at a 
certain threshold Y increases rapidly towards wp. The saturation of cup occurs later 
in Y for increased t (i.e. decreased a;p) as expected from (3.29). Traditional smooth 
diffusion into the infra-red would have led to the opposite behavior, namely the 
higher t case (lower wp) being saturated first. 




5 10 15 20 25 30 

Y 



Figure 7 : The effective intercept c^efr as a function of Y ioT t = 13 and two different values 
of the infra-red cutoff i. 

The study of traditional diffusion is a httle trickier, essentially because for most 
parameter choices there is only a limited domain of traditional diffusion before tunnel- 
ing takes over. Nevertheless one way to approach it is to study the difference between 
the effective exponents in a case with running ccg and one with fixed ctg. For this com- 
parison to be meaningful the fixed ctg must be chosen such that q;sXo(1/2) = 0Js{t)- 
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The observed difference in exponent is then a measure of the difference in the diffu- 
sion properties of the running and fixed-ttg cases: it is sensitive to the terms which 
are expected to be sizeable when a^Y^ ~ 1. 




100 200 300 400 500 600 700 800 900 1000 

Y 

Figure 8: The difference SoJeS between the effective exponents seen with running cxg and 
fixed cts (= 1/101), shown for two values of t. Shown for very extreme values of t and Y 
in order to expose more clearly the main features. 

The difference is shown in Fig. 8, where the dashed line is based on Eq. (3.28), 
which in turn reproduces the results by Kovchegov and Mueller and by Levin [7] 
(with the proper changes for the Xm ^-nd x'L values). The agreement is good, up to 
the critical value Yt beyond which the Pomeron regime takes over. 

In Fig. 8 we have taken extreme values of t and Y, in order to identify the 
running coupling corrections, and we have thus emphasized the difference cup — cUgit) 
also. With realistic parameters and model, the transition between the two regimes 
could be smoother and shifted to relatively higher energies. 

5. Discussion 

The model we have introduced is somewhat intermediate between the full RG- 
improved small- ,t equation and the (Airy) diffusion model with running coupling. 
Compared to the latter, it has the advantage of a correct colhnear behavior, still ad- 
mitting a Schrodinger-type treatment. It appears, therefore, as a useful laboratory 
for testing approximation methods and theoretical prejudices. 

We have already performed several tests. First, we have shown (Sec. 4.1) that 
indeed the a;-expansion of the 7-representation provides a quite good approximation 
of the regular solution, which is thus perturbative, i.e., independent of the boundary 
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conditions in the strong-coupling region. This lends further support to previous 
results [11] on the anomalous dimensions and the hard Pomeron. 

We have further used the present model to analyze the gluon Green's function, 
by testing its factorisation properties, and by providing its explicit expression (Sec. 3) 
which involves the left-regular solution, and is thus dependent on the strong-coupling 
region. In the present model, this dependence occurs through a single scattering 
coefficient S{u), which carries the spectrum and thus the leading Pomeron singularity 
(Eq. (3.2)). It is expected, by the arguments of Rcf. [16], that a similar decomposition 
may be valid in general, perhaps with several scattering coefficients. 

For two-scale processes, it is of particular interest the result of Sees. 3.3 and 
4.2: in an intermediate Y region a perturbative regime exists where the exponent 
u!s{t) is observable, with corrections of type a^Y^, whose size is in agreement with 
Refs. [7]. On the other hand, for large enough energies, Y > t/{uj^ — 0Js{t)), the 
Regge behavior due to the non-perturbative Pomeron takes over with a "tunneling" 
transition to small Vs which may be, in principle, spectacular for sizeable values of 
a;p - uJs{t). 

Several points need further elucidation. First, a realistic evaluation of cjp and 
of the related effective coupling is still needed, in order to put the question of the 
regimes in perspective, and to start making predictions for two-scale processes. In 
particular, the size of uj^ (including possible unitarity corrections) seems to be crucial 
in order to understand how fast and experimentally relevant is the transition to the 
Pomeron regime. 

Furthermore, the relative importance of the perturbative exponents uJs{t) and 
uJc{t) and of uj-p is still to be cleared up in DIS-type processes, where i-evolution 
plays a much more important role. 

For all these questions the present model is likely to provide useful hints and a 
preliminary understanding. 
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